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SUMMARY 


frhe  feasibility  of  using  short-duration  shock -tube  flows  to  accelerate  micron-size 
particles  of  known  mass  and  composition  to  hypervelocities  has  been  studied  both  analyt¬ 
ically  and  experimentally.  Analytic  studies  that  use  the  Stokes  drag  law  corrected  for 
slip  flow  by  the  semiempirical  Cunningham-Millikan  correction  factor  were  made  for 
various  shock  Mach  numbers  to  30,  and  the  character  of  the  acceleration  problem  was 
outlined.  Experimental  studies  were  made  by  accelerating  micron-size  particles  to 
Mach  13  in  a  2.  9-inch  shock  tube.  These  studies  demonstrated  that  short-duration  shock - 
tube  flows  could  be  used  as  a  means  of  accelerati^  micron-size  particles  fo^exposing 
surfaces  to  hypervelocity  impact.  Optical  properties  of  these  exposed  surfaces  were 
quite  reproducibly  affected  by  exposure,  and  an  analytical  model  of  the  surface  degrada¬ 
tion  phenomenon  satisfactorily  related  measured  reduction  of  surface  reflectance  with  the 
total  energy  of  impinging  particles.  The  analysis  suggested  a  means  of  using  accurate 
measurements  in  the  laboratory  to  predict  reflectance  degradation  for  surfaces  under 
actual  meteoroid  dust  exposure  in  space. '  n  k.  \  A  'TP 
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INTRODUCTION 

The  difficulty  of  providing  accurate  estimates  of  the  performance  of  space  vehicle 
components  under  the  unusual  conditions  of  the  space  environment  can  be  reduced  by 
studying  the  effects  of  these  conditions  in  ground  facilities  wherein  the  environment  has 
been  simulated.  Among  the  numerous  facilities  built  for  this  purpose  over  the  last  few 
years  none  has  so  far  included  the  capability  of  providing  simulation  of  the  meteoroid ' 


environment.  The  difficulty  of  providing  such  a  capability  exists  first  because  of  an  in¬ 
complete  knowledge  of  this  environment.  Also,  there  has  been  a  lack  of  satisfactory 
methods  for  accelerating  known  size  particles  to  the  required  relative  speeds  at  impact 
with  the  target. 

Since  the  dangers  to  space  travel  of  this  meteoroid  feature  of  the  space  environment 
are  well  recognized,  many  efforts  have  been  made,  instead,  to  study  the  phenomenon  of 
the  encounter  of  a  single  particle,  moving  at  meteoric  speeds,  with  a  target  representing 
a  vehicle  surface  or  component  wall.  The  purpose  of  such  studies  has  invariably  been 
the  determination  of  the  penetration  laws  for  such  encounters  with  the  express  goal  of 
protecting  various  types  of  fluid  containers  from  catastrophic  puncture.  Data  compiled 
from  many  space  experiments  (refs.  1  and  2)  indicate,  however,  that  the  frequency  of  en¬ 
counter  with  meteoroids  will  increase  very  rapidly  as  the  mass  of  the  meteoroid  particle 
decreases.  For  meteoroid  particle  mass  below  that  which  might  be  considered  a  pene¬ 
tration  hazard,  the  measured  flux  has  increased  so  much  that  the  primary  space  flight 
problem  due  to  these  smaller  meteoroids  becomes  one  of  surface  erosion  and  the  subse¬ 
quent  modification  of  the  optical  properties  of  the  eroded  surfaces.  It  is  clear  that  a  de¬ 
vice  which  could  provide  quantitative  and  reproducible  exposure  of  surfaces  to  high-speed 
micron-size  particles  could  be  useful  in  studying  the  micrometeoroid  erosion  problem, 
and  it  was  with  this  purpose  in  mind  that  a  program  to  develop  such  a  facility  was  initia¬ 
ted. 

The  short -duration  flows  of  high-speed  high-temperature  gas  in  a  shock  tube  have 
been  used  extensively  to  provide  the  environment  necessary  for  the  study  of  the  problems 
of  hypersonic  flight.  The  possibility  of  using  the  aerodynamic  drag  of  such  gas  flows  as  a 
means  of  accelerating  small  particles  of  known  mass  and  composition  to  shocked  gas 
speeds  seemed  worthwhile  to  investigate  because  such  an  investigation  could  be  based  on 
this  strong  aerodynamic  foundation.  It  would  only  involve,  therefore,  the  problem  of 
outlining  the  regions  of  flow  in  the  shock  tube  where  small  particles  could  be  usefully 
accelerated.  An  evaluation  of  the  method  of  exposing  targets  to  particles  accelerated  in 
this  way  could  follow,  hopefully  allowing  the  development  of  a  ground  facility  which  would 
remove  the  need  for  subjecting  every  material  or  surface  of  interest  to  actual  space  ex¬ 
posure. 

Thus,  a  program  was  initiated  at  the  Lewis  Research  Center  to  evaluate  the  use  of 
shock  tubes  for  accelerating  small  particles  of  known  mass  and  composition  to  high 
speeds,  and  further  to  impact  these  high-speed  particles  against  target  surfaces  for  the 
purpose  of  studying  quantitatively  the  effect  of  such  exposure  on  the  optical  properties  of 
these  target  surfaces.  The  studies  made  were  both  experimental  and  analytical.  The  ex¬ 
perimental  studies  were  made  for  various  types  and  size  ranges  of  particles  for  shock 
Mach  numbers  in  air  ranging  from  6. 0  to  13.  0.  The  analytic  studies  were  extended  to 
Mach  30  in  an  effort  to  predict  the  velocities  of  various  size  particles  for  gas  speeds  up 


2 


to  30  000  feet  per  second  in  flow  durations  as  short  as  27  microseconds.  The  develop¬ 
ment  of  the  experimental  methods  and  the  results  of  these  experimental  and  analytical 
studies  are  presented  herein. 


ACCELERATION  OF  PARTICLES 

The  environmental  conditions  encountered  by  a  particle  in  a  shock  tube  can  be  illus¬ 
trated  by  the  wave-particle  diagram  presented  in  figure  1.  The  tube  consists  of  two  cham¬ 
bers,  high  pressure  and  low  pressure,  separated  by  a  diaphragm.  When  the  diaphragm 
opens,  compression  waves,  which  rapidly  steepen  into  a  shock  wave,  move  into  the  low- 
pressure  region,  region  1.  At  the  time  the  diaphragm  opens,  an  expansion  wave  propa¬ 
gates  into  the  high-pressure  driver  gas,  region  4,  reflects  off  the  back  wall  of  the  driver 
section,  and  propagates  into  the  low-pressure  chamber.  The  surface  that  separates 
the  driven  gas  from  the  driver  gas  is  referred  to  as  the  contact  surface.  Across  this 
surface  the  velocity  and  pressure  are  equal,  but  the  temperature  and  density  are  not  in 
general  equal  for  the  two  gases.  (The  temperature  of  the  driven  gas,  air,  at  the  condi¬ 
tions  presented  in  this  report  varies  between  2820°  and  10  875°  K. )  Particles  placed  in 
the  middle  of  the  low-pressure  chamber  at  some  distance  downstream  of  the  diaphragm 
are  accelerated  by  the  driven  gas,  region  2,  and  some  heavier  ones  that  are  not  as  yet  at 
gas  speed  on  arrival  of  the  contact  surface  fall  behind  and  are  accelerated  by  the  cold  ex¬ 
panded  driver  gas,  region  3,  until  the  termination  of  flow,  assumed  here  to  be  the  arrival 
of  the  first  expansion  wave.  It  is  sometimes  desirable  to  have  the  particles  pass  behind 
the  contact  surface  and  continue  to  be  accelerated  by  the  cold  flow  in  region  3  if  the  melt- 


Figure  1.  -  Wave-particle  diagram  for  constant  area  shock  tube. 
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ing  or  sublimation  temperature  of  the  particles  is  less  than  the  hot  gas  temperature  of 
region  2  (since  this  can  produce  an  undesirable  ablating  particle). 

With  the  preceding  flow  picture  in  mind,  an  analysis  was  made  to  determine  the 
feasibility  of  accelerating  particles  with  these  short-duration  shock -tube  flows.  The  ac¬ 
celerating  force  was  assumed  to  be  the  aerodynamic  drag  of  the  high-speed  gas  acting  on 
spheroidal  particles  placed  in  the  stream.  It  was  also  assumed  that  the  time  for  particle 
pickup  at  the  injection  station  was  small,  and  calculations  showed  that  the  Reynolds  num¬ 
bers  for  the  particles  relative  to  the  gas  would  be  low  during  a  large  part  of  the  accel¬ 
erating  process.  These  low  Reynolds  numbers  would  allow  use  (thus  giving  conservative 
drag  estimates)  of  the  well-known  Stokes  solution  for  spheres.  Stokes  drag  has  been  used 
by  many  investigators  in  the  past  and  is  used  here  essentially  to  show  the  character  of 
the  acceleration  problem  in  a  shock  tube.  In  this  case  the  drag  force  is  given  by  the  fol¬ 
lowing  equation  (see  appendix  A  for  definition  of  symbols): 

f^  =  35r?7D(Vg  -  Vp)  (1) 

Now  the  Stokes  drag  law  is  valid  up  to  the  limitations  imposed  by  the  phenomenon  of 
slip.  This  phenomenon  becomes  important  when  the  particle  diameter  is  of  the  same 
order  of  magnitude  as  the  mean  free  path  of  the  molecules  of  the  gas  in  which  the  parti¬ 
cles  are  suspended.  A  correction  for  slip  (ref.  3)  was  made  by  putting  the  semiempirical 
Cunningham-Millikin  correction  factor  k^^  in  the  denominator  of  the  right  side  of  equa¬ 
tion  (1).  For  air  k^^  is  defined  by  the  following  equation: 


^  ^  ^  0. 16xi0~‘^cm  T  ^o  ^2) 
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where  and  T^  are  the  conditions  at  standard  temperature  and  pressure  (STP).  If 
D,  the  diameter  of  the  particle,  becomes  large  enough  to  make  the  second  term  in  equa¬ 
tion  (2)  much  less  than  1,  there  is  no  correction  for  slip  flow.  Since  conditions  in  a 
shock  tube  are  such  that,  at  times,  the  mean  free  path  of  the  accelerating  gas  approaches 
the  particle  diameter,  k^  was  added  to  give  a  conservative  estimate  of  the  acceleration 
process.  Hence,  if  spheroidal  particles  are  assumed,  equation  (1)  can  be  rewritten  as: 


where 
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In  the  present  development,  Vg,  the  gas  velocity,  has  been  assumed  constant.  This 
is,  of  course,  not  strictly  true,  since  the  presence  of  the  particles  (and  the  viscous  ef¬ 
fects  at  the  tube  wall)  act  to  change  the  gas  velocity,  as  well  as  the  shock-wave  velocity, 
as  the  shock  progresses  down  the  tube.  The  phenomena  involved,  however,  do  not  all  act 
to  change  the  gas  velocity  in  the  same  direction  and  have  only  a  small  effect  for  the  quan¬ 
tities  of  material  injected  into  the  size  tube  considered  here.  The  mass  effect  of  the 
particles  is  to  speed  up  the  gas  at  distances  further  from  the  shock  (refs.  4  and  5)  (simi¬ 
lar  to  the  wall  effect  in  ref.  6),  whereas  the  effect  of  the  volume  of  the  particles  is  in  the 
other  direction  (i.  e. ,  to  slow  down  the  gas  behind  the  shock).  In  an  experiment,  more 
material  injected  into  a  smaller  tube  could  produce  a  noticeable  effect,  but  calculations 
made  for  the  particle  masses  considered  herein  indicated  mass  and  volume  loadings  to  be 
sufficiently  low  to  have  a  small  effect  on  the  gas  velocity,  hence  Vg  is  assumed  con¬ 
stant.  Thus,  integration  of  equation  (3)  yields 


V  =  V  + 
P  g 


where  subscript  1  here  refers  to  initial  conditions  on  arrival  of  the  shock  front  at  the 
injection  point,  Xj.  Integration  of  equation  (4)  yields  the  particle  path  equation 


(V„  1  -  V  J 

Xp  =  Xi^Vg(t-t,).  g 


r  -K(t-t )] 

1  -  e  ^ 


(5) 


If  the  particle  is  at  rest  when  the  shock  front  arrives,  and  if  tj^  and  Xj^  are  set  equal 
to  zero,  equations  (4)  and  (5)  reduce  to 


Vp  =  Vg(l  -  e-*^^) 


(6) 
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(7) 


X  =  V  t  -  Xi  (1  - 
P  g  K 

Equations  (6)  and  (7)  were  used  to  calculate  the  velocities  of  various  particles  with  dif¬ 
ferent  values  of  K  (K  varies  for  different  types  and  sizes  of  particles  and  with  the  vis¬ 
cosity  of  the  accelerating  gas)  while  they  remained  in  the  slug  of  shocked  hot  gas.  Parti¬ 
cles  that  are  not  as  yet  up  to  gas  speed  on  the  arrival  of  the  contact  surface  (hot- gas  flow 
duration  was  calculated  by  the  methods  of  ref.  7)  pass  through  the  contact  surface  and  are 
still  accelerated,  but  at  a  slower  rate,  by  the  cold,  expanded  driver  gas. 

After  the  contact  surface  reaches  the  particle,  equations  (4)  and  (5)  were  used  to  cal¬ 
culate  the  trajectory  during  acceleration  by  the  driver  gas.  Since  the  gas  conditions 
{rj  and  k  )  have  changed,  a  new  value  of  K  is  calculated,  X.,  now  becomes  X  ,  V 
equals  the  velocity  of  the  particle  at  the  contact  surface,  and  t^^  here  is  the  time  at 
which  a  particle  meets  the  contact  surface.  The  position  downstream  of  the  injection 
point  at  which  a  particle  meets  the  contact  surface  is  obtained  by  equating  the  distance  the 
contact  surface  has  traveled  from  the  injection  station 

Xc=Vg(t-T)  (8) 

with  the  distance  the  particle  has  traveled  (eq.  (7)) .  It  was  assumed  that  the  injection 
point  of  the  particles  is  far  enough  downstream  of  the  diaphragm  so  that  the  shock -tube 
flow  depicted  in  figure  1  (p.  3)  has  been  established. 


Heating  of  Particle 

The  driven  gas  (air)  behind  a  strong  shock  has  been  elevated  to  high  temperatures. 
For  the  conditions  considered  in  this  report,  real  gas  effects  were  included  to  calculate 
the  temperatures  that  ranged  from  2970°  to  10  875°  K  (ref.  8).  This  hot  gas,  used  to 
accelerate  particles,  can  also  heat  the  particles  to  sufficiently  high  temperatures  to  make 
their  trajectories  visible  on  film.  To  use  a  photographic  method  for  determining  particle 
velocities,  it  was  therefore  necessary  to  determine  the  temperature  history  of  the  parti¬ 
cle.  The  temperature  of  an  accelerating  particle  in  the  hot  gas  at  a  time  t  is  obtained 
by  solving  the  convection  equation: 


-YCo—  =  Ah(T  -  T J 
P  dt  g 


(9) 
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Equation  (9)  is  solved  in  appendix  B  for  the  assumption  that  the  specific  heat  is  con¬ 
stant  over  the  temperature  range  of  interest  and  that  hD/k  =  1  for  a  sphere  with  a 
Reynolds  number  less  than  10  (ref.  9),  The  solution  is 

-t/Cj 

’rp,l  =  Tg+(T^-Tg)e  (10) 

The  film  chosen  for  this  experiment  was  a  high-speed  panchromatic  film  with  a  1250  ASA 
rating. 

The  minimum  energy  density  needed  for  reasonable  exposures  using  this  film  is 
76X10"^  erg  per  square  centimeter  (ref.  10).  To  determine  the  particle  temperature 
necessary  to  attain  this  energy  density,  the  radiation  law  equation  was  solved  in  appen¬ 
dix  B,  with  the  following  result; 


ERVp(t) 

(2.5)2ar|D 


(11) 


The  time  at  which  the  temperature  of  the  particle  in  the  hot  stream  Tp  reaches  the 
temperature  required  to  expose  the  film  Tp  g  was  calculated  from  equations  (10) 
and  (11).  This  time  was  then  put  into  the  position  equation  (7)  to  give  the  mihimum  posi¬ 
tion  downstream  of  the  injection  point  at  which  a  particle  radiates  enough  energy  to  be 
recorded. 

If  the  gas  temperature  T^  is  greater  than  the  sublimation  temperature  of  a  particle, 

any  particle  that  remains  in  the  hot  gas  long  enough  will  reach  its  sublimation  point.  The 

time  t„  at  which  this  takes  place  is  obtained  by  solving  equation  (10)  for  the  time  at 
s 

which  T  1  equals  the  sublimation  temperature  T„: 
p,  i  s 


Equation  (7),  with  t  =  tg  gives  the  position  downstream  of  the  injection  point  at  which 
sublimation  begins. 

Once  a  particle  reaches  the  sublimation  point,  it  continues  to  sublime  until  it  is 
gone,  if  it  does  not  pass  behind  the  contact  surface  into  the  cold  driver  gas.  The  time  at 
which  a  particle  sublimes  to  zero  radius  can  be  obtained  by  solving  the  following  heat- 
balance  equation 
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(13) 


Ah(T„  -  T„)  =  -pEArrr^  — 
g  s  s  dt 


arid  solving  for  t 


t 


s— o 


k(Tg  -  Tg) 


(14) 


Such  a  calculation  was  made  in  the  cases  for  which  it  was  necessary.  The  various  condi¬ 
tions  (initial  pressures  and  gases)  for  which  calculations  were  made  were  chosen  to  be  a 
starting  point  in  picking  optimum  conditions  to  accelerate  particles  in  a  shock  tube. 

Since  a  goal  in  these  experiments  is  the  precise  knowledge  of  particle  size,  the  sublima¬ 
tion  region  should  be  avoided.  The  complication  of  evaluating  the  reduction  in  size  of  the 
particles  due  to  sublimation  was  therefore  not  considered  in  the  present  calculations. 

Such  a  reduction  can  be  included,  if  desired,  by  solving  equation  (13)  for  the  radius  of  the 
particle  and  substituting  r  into  equations  (6)  and  (7)  to  calculate  the  true  trajectory  of 
the  subliming  particle. 


Regimes  for  Particle  Acceleration 

In  most  shock  tubes,  a  limitation  is  imposed  by  the  upper  value  to  which  the  driver 
chamber  can  be  pressurized  or  the  temperature  to  which  the  driver  gas  can  be  raised. 

In  the  analytical  study  to  determine  conditions  in  the  shock  tube  under  which  useful  accel¬ 
eration  of  small  particles  could  be  made,  three  different  driver  gas  conditions  were  con¬ 
sidered,  The  first  two  driver  conditions  require  helium  or  hydrogen  at  room  tempera¬ 
ture,  and  the  third  requires  electrically  heated  helium  (to  give  a  shock  Mach  number 
of  30).  With  driver  chamber  conditions  thus  fixed,  the  maximum  shock  Mach  number 
strongly  depends  on  minimum  pressure  obtainable  in  the  low  pressure  chamber.  How¬ 
ever,  the  actual  duration  of  uniform  flow  of  shocked  gas  also  depends  critically  on  down¬ 
stream  pressures,  decreasing  as  the  pressure  is  decreased.  (Duration  of  uniform  flow 
also  decreases  with  increasing  Mach  number.)  With  the  preceding  considerations  in 
mind,  different  initial  low-pressure  side  pressures  were  chosen  for  the  calculations  to 
give  a  range  of  shock  Mach  numbers  from  6.  5  to  30.  A  variety  of  particles  was  con¬ 
sidered  for  injection  into  and  acceleration  by  the  shocked  gas  flow.  Equations  (4)  to  (7) 
were  used  to  calculate  curves  for  velocity  as  a  function  of  position  downstream  of  injec¬ 
tion.  These  are  presented  as  the  solid  curves  in  figure  2  for  silicon  carbide  (figs.  2(a) 
to  (e)),  tungsten  (figs.  2(f)  and  (g)),  and  graphite  (figs.  2(h)  and  (i))  particles  to  show  the 
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character  of  the  problem  of  accelerating  such  particles  in  these  short-duration  flows. 
Shown  as  dot-dashed  curves  superimposed  over  these  curves  are  the  boundary  curves  for 
the  points  beyond  which  the  various  particles  are  visible  (as  calculated  from  eqs.  (7). 

(10),  and  (11)).  The  dashed  curve  superimposed  in  each  of  these  figures  is  the  locus  of 
positions  at  which  the  various  particles  meet  the  contact  surface  (eq.  (8)).  The  locations 
at  which  the  particles  first  begin  to  sublime,  are  joined  by  a  dash-dash-dot-dot  curve  in¬ 
dicating  the  boundary  beyond  which  the  particle  size  would  vary  from  its  original  value 
(eqs.  (7)  and  (12)). 

For  the  silicon  carbide  particles  ranging  from  2  to  20  microns  and  accelerated  by  gas 
behind  a  shock  of  M  =  6.  5  (fig.  2(a)),  only  the  smaller  particles  (<10|U,)  become  visible 
and  are  accelerated  to  near  gas  speed  (horizontal  dashed  line  in  the  figure)  before  reach¬ 
ing  the  cutoff  point  (in  these  figures,  this  point  is  10  ft  downstream  of  injection).  The  14- 
and  20-micron  particles  pass  behind  the  contact  surface  before  becoming  visible,  and 
hence  never  become  visible  during  acceleration  by  the  cold  expanded  driver  gas.  The  uni¬ 
form  shocked-gas  flow  duration  for  this  case  was  375  microseconds.  A  2-micron  particle 
starts  to  sublime  at  0.  557  foot,  and  larger  particles  would  start  to  sublime  further  down¬ 
stream  (5-/i  particle  at  4. 08  ft).  Thus,  in  this  figure,  the  sublimation  curve  collapses 
onto  the  particle  velocity  curves  near  gas  velocity  and  is  not  shown. 

In  figure  2(b),  the  shock  Mach  number  has  been  increased  to  8.  38  by  reducing  the 
downstream  pressure.  This,  as  mentioned  previously,  causes  the  flow  duration  to  de¬ 
crease,  in  this  case  to  142  microseconds.  The  effect  immediately  apparent  is  the  re¬ 
duction  in  the  range  of  size  of  particles  (<5p.),  which  can  be  accelerated  to  near  gas  speed 
in  the  cutoff  distance. 

In  figure  2(c),  the  initial  pressure  downstream  is  reduced  to  1  millimeter  of  mer¬ 
cury,  thereby  reducing  the  accelerating  flow  duration  to  98  microseconds  while  only  in¬ 
creasing  the  shock  Mach  number  slightly.  The  reduction  in  flow  duration  reduces  further 
the  size  of  particle  that  can  be  accelerated  to  near  gas  speed  by  the  cutoff  distance.  Here 
there  is  no  advantage  to  be  gained  in  decreasing  the  downstream  pressure,  since  shock- 
tube  pressure-ratio  limitations  on  Mach  number  are  being  reached  and  the  only  effect  is 
reduction  in  flow  duration. 

If  a  driver  gas  of  lower  molecular  weight  is  used,  with  no  change  in  downstream 
pressure,  a  higher  shock  Mach  number  is  attainable,  hi  figure  2(d)  the  driver  gas  is  hy¬ 
drogen,  and  a  gas  speed  of  13  000  feet  per  second  is  obtained.  The  higher  accelerating 
gas  temperatures  associated  with  the  increased  Mach  number  cause  the  gas  to  have  higher 
viscosity  and  hence  provide  higher  drag.  However,  earlier  sublimation  due  to  the  higher 
temperatures  limits  the  minimum  size  particle  that  can  be  used,  and  the  shorter  flow 
duration  (43/isec)  at  the  higher  Mach  number  limits  the  maximum  particle  size.  Clearly, 
the  range  of  size  of  particles,  which  can  be  usefully  accelerated  by  the  shocked  gas  is 
diminishing,  and  there  appears  to  be  a  near  upper  Mach  number  limit  for  accelerating 
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Particle  velocity,  Vp,  ft/sec 


(a)  Silicon  carbide  particles  inserted  20  feet  from  diaphragm.  Driver  gas,  helium  at  1000  pounds  per 
square  inch  pressure;  driven  gas,  air;  initial  driven  gas  pressure,  5. 7  millimeters  of  mercury; 
Mach  number,  6. 5;  gas  velocity,  6050  feet  per  second;  gas  temperature  in  region  2,  3000°  K;  flow 
duration,  375  microseconds. 


(b)  Silicon  carbide  particles  inserted  20  feet  from  diaphragm.  Driver  gas,  helium  at  1000  pounds  per 
square  inch  pressure;  driven  gas,  air;  initial  driven  gas  pressure,  2.0  millimeters  of  mercury;  Mach 
number,  8. 38;  gas  velocity,  8250  feet  per  second;  gas  temperature  in  region  2,  2970°  K;  flow  duration, 
142  microseconds. 

Figure  2.  -  Particle  velocity  against  distance  downstream  from  injection  point 
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Particle  velocity,  Vp,  ft/sec 
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(c)  Silicon  carbide  particles  inserted  20  feet  from  diaphragm.  Driver  gas,  helium  at  1000  pounds  per 
square  inch  pressure;  driven  gas,  air;  initial  driven  gas  pressure,  1.0  millimeter  of  mercury;  Mach 
number,  8.7;  gas  velocity,  8410  feet  per  second;  gas  temperature,  in  region  2,  3000°  K;  flow  duration 
98  microseconds.  \  . 
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Distance  from  injection  point  of  particle,  Xp,  ft 

(d)  Silicon  carbide  particles  inserted  20  feet  from  diaphragm.  Driver  gas,  hydrogen  at  1500  pounds  per 
square  inch  pressure;  driven  gas,  air;  initial  driven  gas  pressure,  1.0  millimeter  of  mercury;  Mach 
number,  13. 3;  gas  velocity,  13  100  feet  per  second;  gas  temperature  in  region  2,  4890°  K;  flow  dura¬ 
tion,  43  microseconds. 


Figure  2.  -  Continued. 


Particle  velocity,  ft/sec 


(e)  Silicon  carbide  particles  Inserted  60  feet  from  diaphragm.  Driver  gas,  helium  at  dOOO  pounds  per  square  inch 
pressure  electrically  heated  to  16  900°  K;  driven  gas,  air;  initial  driven  gas  pressure,  1. 0  millimeter  of  mercury; 
Mach  number,  30;  gas  velocity,  30  250  feet  per  second;  gas  temperature  in  region  2,  10  875°  K;  flow  duration, 

31  microseconds. 


If)  Tungsten  particles  inserted  20  feet  from  diaphragm.  Driver  gas,  helium  at  1000  pounds  per  square  inch  pres¬ 
sure;  driven  gas,  air;  initial  driven  gas  pressure,  6. 7  millimeters  of  mercury;  Mach  number,  6. 5;  gas  velocity, 
6050  feet  per  second;  gas  temperature  in  region  2,  3000°  K;  flow  duration,  375  microseconds. 


Figure  2.  -  Continued. 
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Particle  velocity,  Vp,  ft/sec 


(g)  Tungsten  particles  inserted  30  feet  from  diaphragm.  Driver  gas,  helium  at  1000  pounds  per  square  inch  pres¬ 
sure;  driven  gas,  air;  initial  driven  gas  pressure,  6.7  millimeters  of  mercury;  Mach  number,  6.5;  gas  velocity, 
6050  feet  per  second;  gas  temperature  in  region  2,  3000°  K;  flow  duration,  507  microseconds. 


Distance  from  injection  point  of  particle,  Xp,  ft 

(h)  Graphite  particles  inserted  20  feet  from  diaphragm.  Driver  gas,  hydrogen  at  1500  pounds  per  square  inch  pres¬ 
sure;  driven  gas,  air;  initial  driven  gas  pressure,  1. 0  millimeter  of  mercury;  Mach  number,  13.3;  gas  velocity, 
13  100  feet  per  second;  gas  temperature  in  region  2,  4890°  K;  flow  duration,  43  microseconds. 

Figure  Z  -  Continued. 
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particles  of  a  given  material  with  shock -tube  flows.  However,  in  attempting  to  obtain 
even  higher  shock  Mach  numbers,  shock-tube  pressure- ratio  limitations  are  overcome  by 
heating  the  driver  gas  (ref,  11),  This  heating  (to  16  900°  K  as  in  fig.  2(e))  not  only  in¬ 
creases  the  Mach  number  of  the  shock,  but  also  increases  the  effectiveness  of  the  ex¬ 
panded  driver  gas  behind  the  contact  surface.  As  shown  in  figure  2(e),  this  driver  gas 
(which  is  at  2050°  K  after  expansion)  quite  effectively  accelerates  a  range  of  particle 
sizes  after  the  particles  have  passed  through  the  contact  surface.  Particles  5,  6,  and  7 
microns  in  diameter  remain  visible  and  are  accelerated  to  near  gas  speeds  before  reach¬ 
ing  the  cutoff  point  at  10  feet.  This  increase  in  effectiveness  is  a  result  of  the  increased 
viscosity  of  the  heated  driver  gas  even  after  expansion.  Thus,  the  apparent  limitation  on 
accelerating  particles  to  higher  speeds  in  shock -tube  flows  may  be  overcome  by  heating 
the  driver  gas. 

Further  calculations  were  made  for  tungsten  particles,  to  see  whether  material 
properties  could  be  used  for  avoiding  the  sublimation  limitations  previously  obtained. 
Figure  2(f)  shows  that,  because  of  the  higher  sublimation  temperature  of  tungsten,  no 
sublimation  occurs,  but  the  larger  density  imposes  serious  limitations  on  the  maximum 
size  of  particle  that  can  be  accelerated  to  near  gas  speed.  Injecting  particles  further 
downstream  increases  the  hot-gas  flow  duration,  and  in  figure  2(g),  the  curves  for  tung¬ 
sten  particles  injected  30  feet  downstream  of  the  diaphragm  show  that  the  maximum  size 
of  particles  usefully  accelerated  does  increase.  However,  inserting  particles  further 
downstream  in  a  given  length  tube  automatically  reduces  the  length  of  tube  remaining  af¬ 
ter  injection.  Hence,  this  mechanism  for  increasing  acceleration  time  is  not  necessarily 
advantageous. 

To  give  maximum  range  in  size  of  particles  that  can  be  accelerated  to  near  gas 
speed,  the  particle  material  selected  should  have  a  high  sublimation  point  and  a  low  den¬ 
sity.  These  criteria  are  borne  out  by  the  calculations  presented  in  figures  2(a)  to  (g). 

To  make  the  particles  visible  for  speed  measurements,  however,  specific  heat  is  an  addi¬ 
tional  characteristic  that  must  be  considered. 

In  figure  2(h),  curves  for  carbon  (graphite)  particles  show  that  material  with  a  rela¬ 
tively  high  sublimation  temperature  and  a  relatively  low  density  is  still  not  ideal  because 
of  relatively  low  specific  heat.  Although  such  particles  become  visible  rapidly,  they  be¬ 
come  invisible  rapidly  after  passing  behind  the  contact  surface.  Thus  choice  of  materials 
eases  but  does  not,  in  general,  remove  the  problems  encountered. 

Heating  of  the  driver  gas  appears  then  to  be  the  most  effective  means  of  removing 
limitations  on  microparticle  acceleration  in  shock-tube  flows.  This  conclusion  is  further 
ther  strengthened  by  the  calculations  for  carbon  particles  shown  in  figure  2(i),  in  which 
the  heated  helium  driver  gas  successfully  accelerates  6-  and  7-micron-diameter  graphite 
particles  to  30  250  feet  per  second  while  keeping  the  particles  visible  and  avoiding  the 
problems  of  sublimation. 
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(i)  Graphite  particles  inserted  40  feet  from  diaphragm.  Driver  gas,  helium  at  4000  pounds  per  square  inch  pres¬ 
sure  electrically  heated  to  16  900°  K;  driven  gas,  air;  initial  driven  gas  pressure,  1. 0  millimeter  of  mercury,- 
Mach  number,  30;  gas  velocity,  30  250  feet  per  second;  gas  temperature  in  region  2,  10  875°  K;  flow  duration, 

27  microseconds. 

Figure  2.  -  Concluded. 

The  preceding  calculations  were  made  by  using  the  Stokes  drag  law,  although  this 
is  conservative,  to  show  the  character  of  the  problems  involved  when  accelerating  micro¬ 
particles  in  shock-tube  flows.  Actual  drag  coefficients  may  be  as  much  as  an  order  of 
magnitude  larger  than  the  Stokes  value  (refs.  12  and  13).  This  would  not  change  any  con¬ 
clusions  drawn  from  the  material  presented  in  figure  2.  It  would  only  have  the  effect  of 
making  the  curves  shown  valid  for  slightly  larger  particles. 

LABORATORY  EXPOSURE  OF  REFLECTIVE  SURFACES 

The  short-duration  flows  in  a  2.  9-inch-diameter  shock  tube  (shown  in  fig.  3(a))  with 
a  driven  section  35  feet  long  (fig.  3(b))  were  used  to  accelerate  micron-size  particles  in 
the  experimental  portion  of  the  program.  Helium  at  1000  poimds  per  square  inch  was 
used  as  the  driver  gas,  and  air  was  the  driven  gas  in  the  low-pressure  section  into  which 
the  shock  was  fired.  A  copper  diaphragm  separating  the  chambers  was  ruptured  by  an 
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(a)  Overall  view. 


Helium 


(b)  Schematic  drawing  indicating  position  of  shock  velocity  pickups  and  camera  iocation. 
Figure  3.  -  2.9-lnch-diameter  shock  tube. 


16 


arrowhead  plunger.  This  system  was  very  successful  in  initiating  shocked  flows  while 
preventing  any  diaphragm  particles  from  tearing  loose  and  interfering  with  the  experi¬ 
ment. 

The  driven  end  of  the  2,  9-inch  shock  tube  consisted  of  two  tube  sections  in  series, 

20  feet  long  and  9.  6  feet  long,  respectively.  These  were  followed  by  a  transition  section 
10  inches  long  and  a  rectangular  test  section  4.  46  feet  long  which  emptied  into  an  18-inch- 
diameter  dump  tank  8.  6  feet  long.  Shock-speed  measurements  were  made  with  thin-film 
platinum  gages  on  the  wall  of  the  driven  tube  just  upstream  of  the  transition  section. 
Particle  velocity  measurements  were  made  by  streak  film  photography  of  the  glowing 
particles  in  the  test  section.  The  particle  collection  station  was  just  downstream  of  the 
test  section  in  the  dump  tank.  Test  surfaces  placed  at  the  collection  station  were  exposed 
to  impaction  by  a  known  mass  of  particles  arriving  at  measured  speeds.  For  actual  ex¬ 
posures  of  surfaces,  only  particles  that  would  be  accelerated  to  gas  speeds  were  used  in 
the  experiment.  This  exposure  was  correlated  with  the  change  of  test  surface  properties 
as  measured  by  optical  instruments.  The  optical  measurements  made  were  of  test  sur¬ 
face  reflectance  before  and  after  exposure  and  involved  the  use  of  a  spectral  reflectom- 
eter  in  conjunction  with  a  heated  cavity  (hohlraum)  blackbody  radiator. 


Particle  Injection 

Particles  to  be  injected  into  the  stream  were  placed  on  the  horizontal  surface  of  a 
sharp-edged  thin  plate  located  at  the  injection  station  in  the  center  of  the  tube  cross  sec¬ 
tion.  It  was  from  this  surface  that  the  particles  were  picked  up  by  the  high-speed  flows 
following  passage  of  the  shock.  The  data  presented  herein  were  obtained  for  silicon  car¬ 
bide  particles  injected  at  a  point  20  feet  downstream  of  the  diaphragm.  The  injection  sys¬ 
tem  required  great  care  but  with  this  care  quite  repeatable  particle  measurements  of 
collector  weights  and  distribution  of  particles  were  possible. 


Velocity  Measurements 

Velocity  measurements  of  the  gas-heated,  radiating  cloud  of  particles  were  made 

by  recording  particle  trajectories  with  a  streak  camera  on  1250  ASA  film.  The  camera, 

which  had  a  depth  of  field  of  1/8  inch,  was  focused  on  the  center  vertical  plane  in  the  test 

section  of  the  shock  tube  (see  fig.  3(b)).  This  plane  is  parallel  to  the  longitudinal  axis  of 

the  tube  and  therefore,  of  course,  also  parallel  to  the  flow  direction.  The  vertical  sides 

of  the  rectangular  test  section  of  the  shock  tube  are  provided  with  windows  that  allow  ob- 

5  5 

servation  of  the  flow  of  particles.  A  slit  2g  inches  long  and  g  inch  high  was  mounted  on 
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Figure  4.  -  Streaks  made  by  radiation  given  off  by  2-  to  14- 
micron-diameter  silicon  carbide  particles  at  Mach  number  6.5 
taken  4.45  feet  downstream  of  injection  station. 


the  front  of  one  of  the  glass  windows.  A  typ¬ 
ical  photograph  is  shown  in  figure  4.  In  this 
case,  the  streaks  were  made  by  radiation 
from  silicon  carbide  particles  2  to  14  mi¬ 
crons  in  diameter.  The  shock  Mach  number 
was  6.  5,  and  the  initial  pressure  Pj  was 
6.  7  millimeters  of  mercury.  The  velocity  is 
determined  from  the  angle  a  streak  on  the 
film  makes  with  the  horizontal.  The  speed 
of  the  film  (the  time  axis)  is  accurately 
known  (to  a  small  fraction  of  a  percent)  from 
precise  measurements  of  the  rotational  speed 
of  the  mirror,  which  beams  the  picture  to 
the  film  at  the  focal  cylinder  in  the  camera 
(mirror  rotational  speeds  are  240  to  500 
rps).  Thus,  particle  speeds  can  be  mea¬ 
sured  with  good  accuracy.  To  make  certain 
that  the  streaks  used  for  velocity  measure¬ 
ments  were  made  by  particles  and  not  by 
glowing  gas,  runs  were  made  for  which  no 
particles  were  injected  into  the  stream  and 
the  film  exposed  in  the  usual  way.  Invari¬ 
ably,  after  such  runs  the  film  appeared  un¬ 
exposed,  indicating  that  particles  were  neces¬ 
sary  to  obtain  the  streaks  on  the  film,  and 
that,  therefore,  the  streaks  were  indeed 
made  by  radiating  particles.  Velocity  mea¬ 
surements  were  also  successfully  made  for 
aluminum  oxide,  nickel  and  carbon,  as  well 
as  for  glass  microballoon  particles,  although 
these  were  not  used  for  surface  exposures  in 
the  present  study. 

No  attempt  was  made  to  measure  the 
velocity  of  a  single  particle  by  using  the 
streak  camera  technique,  since  the  develop¬ 
ment  of  the  technology  required  to  inject  and 
locate  a  single  particle  would  multiply  the 
difficulties  of  the  experiment  many  times. 
Cases  for  actual  exposures  of  surfaces  were 
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limited  to  those  for  which  all  particles  involved 
would  be  near  gas  speed.  Since  the  original  pur¬ 
pose  of  using  the  shock  tube  was  to  accelerate 
micron-size  particles  for  surface  erosion  studies, 
and  such  erosion  requires  exposure  to  impact  by 
many  particles,  experiments  with  a  single  parti¬ 
cle  were  not  required. 

Distribution  of  Particles  at  Impact 

To  determine  the  distribution  of  the  particles 
over  the  plane  at  the  end  of  the  shock  tube,  mea¬ 
surements  were  made  by  collecting  silicon  car¬ 
bide  particles  on  various  diameter  disk  collectors, 

shown  in  figure  5.  The  collectors  range  in  diam- 
1  13 

eter  from  -g  to  inches  and  were  coated  with 

vacuum  grease  to  ensure  capture  of  the  accelerated  particles.  The  largest  diameter  disk 
13 

(3-jg  in. )  was  held  in  the  tube  by  sting  mounts  and  placed  in  the  dump  tank  4  inches  from 

the  end  of  the  test  section  (see  fig.  3(b),  p.  16).  It  was  presumed,  and  later  verified, 

that  all  the  particles  actually  accelerated  would  be  collected  by  this  disk.  The  smaller 

size  disks  were  placed  in  the  shock  tube  by  mounting  them  on  a  3^ -inch-diameter  plate 

that  was  sting  mounted  in  the  shock  tube  in  the  same  manner  and  at  the  same  station  as 
13 

the  3-jg -inch  collector  disk.  The  coated  disk  collectors  were  weighed  on  an  analytical 
balance  before  and  after  bombardment  by  particles.  The  difference  in  weight  gives  the 
mass  of  particles  collected  by  each  disk,  and  over  the  range  of  collector  disk  size,  the 
radial  density  distribution.  When  no  particles  were  injected,  as  accurately  as  could  be 
measured  (50  jug)  the  weight  of  the  collectors  was  unchanged.  The  mass  of  the  silicon 
carbide  particles  collected,  normalized  to  the  mass  of  particles  collected  on  the  largest 
size  disk  (i.  e. ,  to  the  total  amount  actually  accelerated  down  the  tube)  is  plotted  against 
the  diameter  of  the  disk  collectors  in  figure  6.  It  can  be  seen  that  the  radial  distribution 
of  particles  varied  slightly  with  the  Mach  number  for  those  shown  and  was  essentially 
independent  of  the  total  initial  mass  of  particles  placed  in  the  shock  tube  over  the  range 
of  mass  increments  used  in  this  experiment.  For  any  intermediate  size  disk,  the  mass 
of  silicon  carbide  particles  falling  on  the  disk  could  be  quite  satisfactorily  repeated  to 
within  a  few  percent.  Collections  attempted  for  other  projectile  materials  and  particle- 
size  range  indicated  that  such  accurate  collections  as  with  silicon  carbide  (of  6-fj.  average 
diameter)  could  not  always  be  accomplished  by  the  technique  described.  Some  of  the 
other  types  of  particles  successfully  accelerated  tended  to  agglomerate  before  insertion 
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Figure  5.  -  Disk  coliectors  of  various  sizes  (1/2  to 
13 

3JJ  in.  in  diam)  used  to  make  integrated  density 
measurements. 
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Disk  diameter,  in. 

Figure  6.  -  Mass  of  particles  collected  on  disks  of  various  diameters 
normalized  to  total  mass  collected  as  function  of  diameter  of  disk 
collector. 

in  the  shock  tube  as  a  result  of  collected  moisture  or  static  charge.  These  agglomerates 

removed  large  amounts  of  the  vacuum  grease  coating  the  disk  collectors  and  thus  made 

weight  measurements  with  this  type  of  collector  imfeasible.  For  quantitative  exposures, 

an  arbitrary  minimum  of  75  percent  was  set  for  the  percentage  of  the  total  mass  injected 

13 

to  be  collected  on  the  largest  (3jg  in. )  disk. 

This  minimum  figure  was  set  in  order  to  include  the  possibility  that  some  of  the  mea¬ 
surement  decrement  be  a  result  of  loss  of  collector  mass  and  in  this  way  limit  the  maxi¬ 
mum  error  that  might  arise.  Since  collection  measurements  were  repeatable  to  ±5  percent 
and  more  than  80  percent  of  the  total  mass  injected  was  usually  collected,  this  total 
amount  collected  was  considered  to  be  the  amount  actually  accelerated.  No  total  collection 
weight  measurements  ever  exceeded  the  weight  of  projectile  material  placed  in  the  tube. 
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K  this  minimum  was  not  met  for  a  given  projectile  material  and  size,  no  quantitative 
exposures  were  made.  The  only  projectile  particles  used  for  quantitative  surface  expo¬ 
sure  in  the  present  report  (SiC,  2  to  14  fi)  met  this  requirement.  More  than  75  percent 
of  those  placed  at  the  injection  station  were  collected  on  the  largest  disk.  The  normali¬ 
zation  in  figure  6  is  based  on  the  assumption  that  the  total  mass  actually  accelerated  will 

13 

be  collected  on  the  Sjg -inch-diameter  disk.  The  assumption  appears  to  be  true  since  the 
3 -inch-diameter  disk  collects  almost  the  same  amount  (~95  percent)  as  does  the  larger 
disk.  The  abrupt  leveling  off  of  the  curve  at  3  inches  in  diameter  indicates  that  further 
increase  in  collection  disk  diameter  would  not  increase  the  amount  collected. 


Optical  Properties  of  Impacted  Surfaces 

If  a  polished  metal  surface  is  hit  by  a  microparticle  at  hypervelocity,  a  hemispherical 
crater  is  left  in  the  surface.  The  volume  of  the  crater  depends  on  the  target  material  and 
has  been  shown  to  vary  directly  with  the  kinetic  energy  of  the  projectile  (in  the  range  of 
speeds  considered  here,  ref.  14).  The  area  of  the  surface  damaged  depends,  of  course, 
on  the  crater  volume.  The  resulting  total  reflectance  of  a  surface  can  be  described  in 
terms  of  a  replacement,  over  the  area  of  damage,  of  the  original  reflectance  with  a 
new  reflectance  p^.  This  can  be  stated  mathematically  as 

Pf{AD(e)}  =  ^ 
o 

With  this  model,  it  is  clear  that  further  damage  to  the  surface  can  occur  only  on  the  re¬ 
maining  undamaged  area.  From  this  it  follows  that 

^  -  Acte)]  (16) 

Solving  equation  (16)  for  Ap(e)  gives 

=  ^o  (^  -  ® 

Substituting  for  Aj^(e)  from  equation  (17)  into  equation  (15)  yields 


Ao  -  Aj,(e)  (1  - 


(15) 
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(18) 


Pf(e)  =  Pi 


Equation  (18)  describes  the  reduction  in  reflectance  of  any  exposed  area  that  is 
struck  by  hypervelocity  particles,  the  sum  of  whose  kinetic  energy  is  e.  The  constant 
Kj,  which  is  the  fraction  of  fresh  area  damaged  per  unit  of  exposure  energy,  can  be  ob¬ 
tained  from  equation  (16)  for  the  case  of  zero  exposure  (Aj^  =  0)  as 


(19) 


For  the  first  increment  of  exposure  (i.  e. ,  the  first  hit),  if  the  crater  volume  is  pro¬ 
portional  to  particle  kinetic  energy. 


dA 


D 


de 


-17-2 

Vp 


(20) 


2 

where  the  denominator  m  V^/2  is  the  single  particle  kinetic  energy,  the  numerator  is 

P  P 

the  surface  area  of  the  first  crater,  and  is  the  cratering  energy  density,  that  is, 
the  projectile  kinetic  energy  required  per  unit  volume  of  the  crater  formed  in  the  target. 
The  cratering  energy  density  E^^  has  been  evaluated  for  many  materials  and  is  available 
in  the  literature  (refs.  14  to  16).  Solving  equations  (19)  and  (20)  for  results  in 


K 


1 


(21) 


It  is  clear  from  equation  (21)  that  the  assumption  made  earlier,  that  is  constant  for 
a  given  particle  kinetic  energy,  is  essentially  the  assumption  that  E^^  is  constant. 
Finally,  equation  (18)  becomes 
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Reflectance  as  a  function  of  exposure  can  be  plotted  directly  from  equation  (22)  for  a 
given  target  material,  projectile,  and  projectile  kinetic  energy.  Values  for  E  were 
obtained  from  reference  16,  while  was  the  reflectance  determined  by  exposing  a  tar¬ 
get  to  impaction  by  microparticles  until  no  further  change  in  reflectance  could  be  mea¬ 
sured.  Equation  (22)  was  used  for  comparison  with  the  reflectance  measurements  de¬ 
scribed  in  the  next  section. 


Measurement  of  Optical  Properties 


Test  surfaces  (for  cases  described  herein,  generally  polished  metals)  of  disks  15/16 
inch  in  diameter  and  1/ 64  to  1/16  inch  thick  were  exposed  to  impaction  by  known  amounts 
of  particles  of  known  size  and  composition  at  measured  speeds.  With  exposure  thus 
known,  it  was  desirable  to  determine  whether  the  effects  of  this  exposure  on  surface  op¬ 
tical  properties  would  be  reproducible  at  any  given  exposure.  Infrared  reflectance  was 
chosen  as  an  operational  parameter,  and  spectral  reflectance  measurements  of  exposed 
disks  were  made  from  1.  5  to  15.  5  microns  with  an  infrared  spectrometer  coupled  with  a 
heated-cavity  blackbody.  In  such  a  system,  the  spectrometer  compares  the  radiation  at  a 
given  wavelength  from  a  600°  C  blackbody  cavity  with  the  total  radiation  reflected  from  a 
water-cooled  sample  within  the  cavity.  The  ratio  of  the  amount  reflected  to  the  blackbody 
value  is  taken  as  the  reflectance  at  that  wavelength  and  is  designated  as  ^  To  de¬ 

termine  a  total  (wavelength  independent)  reflectance,  this  spectral  value  is  weighted  for 
the  radiation  from  a  blackbody  at  420°  K.  Total  final  reflectance  is  defined  here  as 


L 


Pf  = 


fh 


(23) 


for  Xj  =  1.  5  microns  and  Xg  =  ^  microns.  It  is  this  average  value  that  characterizes 

the  optical  condition  of  the  surface  and  is  used  to  determine  the  effect  of  exposure. 

^These  data  were  obtained  by  E.  Anagnostou  of  Lewis. 
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Figure  7.  -  Sample  aluminum  disk  15/16  inch  in  diameter  before  and  after  ex¬ 
posure  to  2-  to  14-micron  silicon  carbide  particle  flows. 


Figure  8.  -  Photomicrograph  of  exposed  surface  shown  in  figure  7.  X500. 
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Figure  9.  -  Spectral  reflectance  pj^.g  as  function  of  wavelength  tor  polished  hard  aluminum  (2024)  disk  exposed  to  1.21  joules  of  6-micron 
diameter  siiicon  carbide  particles. 


o  Experimental  data;  - 


/•15.5 
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where  ^(W  =  420°  K  blackbody 
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Figure  10.  -  Plot  of  reflectance  ratio  as  function  of  kinetic  energy  of  6-micron-diameter  silicon  carbide  particies 
falling  on  an  aluminum  disk  15/16  inch  in  diameter.  Total  effective  reflectance  for  area  damaged  by  crater,  poo, 
0. 34;  total  initial  reflectance  of  area  Aq  before  exposure,  Pj,  0. 991;  disk  area,  Aq,  wr^  =  7r(15/32  in.)^;  cratering 
energy  density.  Err,  1. 65x10^0  erg  per  cubic  centimeter;  particle  mass,  m.,  3, 62x10" gram;  particle  velocity, 
V„,  8200  feet  per  sUd.  ^ 


Degradation  of  Reflectance  with  Exposure 

A  photograph  of  a  typical  aluminum  disk  before  and  after  exposure  to  microparticle 
impact  is  shown  in  figure  7.  The  exposed  surface  magnified  500  times  is  shown  in  fig¬ 
ure  8,  where  the  hemispherical  character  of  the  craters  is  evident. 

The  values  of  as  a  function  of  wavelength  for  a  15/ 16-inch-diameter  aluminum 

disk  before  and  after  exposure  to  1.  21  joules  of  silicon  carbide  particles  is  presented  in 
figure  9.  Note  that  the  aluminum  before  exposure  has,  within  the  accuracy  of  the  mea¬ 
surement,  a  reflectance  of  approximately  1.  Dover  the  range  1.  5  to  15.  5  microns.  Averag¬ 
ing  this  with  equation  (23)  gives  a  of  0. 991,  and  this  is  shown  in  the  figure.  After 
exposure,  the  reduction  in  reflectance  is  not  the  same  at  all  wavelengths,  and  the  weight¬ 
ing  by  means  of  equation  (23)  gives  a  value  for  of  0.  793.  Note  that  for  the  shorter 
wavelengths  in  the  infrared  region,  which  are  of  greater  interest  for  solar  radiation,  the 
effect  of  degradation  of  reflectance  is  the  largest. 

The  measured  ratio  of  p^/p-  for  various  exposures  is  presented  in  figure  10,  where 
the  point  for  1.  21  joules  is  this  ratio  for  the  data  presented  in  figure  9.  The  solid  line  is 
a  plot  of  equation  (22)  with  p^  as  measvired  in  the  present  experiment  and  from 
reference  16.  The  excellent  fit  may  be  fortuitous  and  is  not  necessarily  indicative  of  the 
accuracy  to  be  expected  in  this  work.  It  does  indicate,  however,  that  with  reasonable 
values  for  the  material  constants,  a  good  estimate  can  be  made  of  the  surface  degradation 
as  a  result  of  hypervelocity  exposure  to  microparticles  in  the  laboratory. 

CONCLUDING  REMARKS 

An  analysis  is  presented  to  establish  the  feasibility  of  accelerating  microparticles  to 
hypervelocities  in  shock-tube  flows.  Conservative  assumptions  with  respect  to  drag  were 
made  because  they  were  simple,  and  allowed  the  greatest  clarity  in  the  presentation  of 
the  character  of  the  problem.  Regimes  for  useful  acceleration  of  microparticles  were 
found  indicating  that  these  short-duration  flows  can  accelerate  particles  over  a  useful 
range  of  sizes  to  speeds  of  interest.  Experiments  showed  that  sufficiently  accurate 
measurements  of  speed  and  quantity  of  impinging  material  can  be  made  for  characterizing 
the  exposure  of  target  materials.  At  least  for  polished  metal  surfaces,  optical  properties 
of  the  surfaces  correlated  with  the  exposure.  An  analytical  determination  of  the  variation 
in  surface  reflectance  of  these  targets  with  exposure  compares  quite  adequately  with  ex¬ 
perimentally  determined  values.  The  analysis  suggests  a  means  of  extending  measure¬ 
ments  of  reflectance  degradation  in  the  laboratory  to  predictions  for  surfaces  under  actual 
meteoroid  dust  exposure  in  space.  Such  an  estimate  was  made  in  reference  17. 

Lewis  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Cleveland,  Ohio,  September  29,  1965. 


26 


APPENDIX  A 


SYMBOLS 


A 

surface  area  of  particle 

k 

m 

Cunningham- Millikan  correc¬ 

damaged  area  after  exposure 

to  e 

tion  factor, 

^  1  0.  lexlO"*^  cm  T  ^0 

Af 

area  of  film  exposed  by  radia¬ 
ting  particle,  A/ 

D  P 

disk  area 

L 

length  in  object  plane 

b 

3ffi?D/k^ 

Lf 

M 

length  in  focal  plane 

Mach  number 

Cp 

specific  heat  of  particle 

s 

2r2Cpp/3k 

“p 

particle  mass 

D 

particle  diameter 

P 

pressure 

minimum  energy  density 
needed  to  expose  high-speed 

Pi 

initial  downstream  pres¬ 

E 

sure,  mm  Hg 

panchromatic  film  (1250  ASA) 

P4 

initial  driver  pressure 

E 

cr 

cratering  energy  density 

R 

distance  from  particle  to  lens 

F 

demagnification  factor,  L/Lj 

of  streak  camera 

% 

drag  force 

r 

particle  radius 

Hs 

heat  of  sublimation 

radius  of  camera  lens 

h 

heat -transfer  coefficient 

T 

temperature 

spectral  energy  distribution 

gas  temperature 

of  420°  K  blackbody 

^P,l 

temperature  of  particle  in  hot 

K 

b/m  =  ISrj/k  pD^ 

gas 

^1 

fraction  of  fresh  area  damaged 
per  unit  of  exposure  energy. 

Vs 

temperature  required  for 
particle  traveling  at  V  to 
radiate  minimum  energy 
needed  to  expose  film 

AoKr/  V  2  J 

Tr 

room  temperature 

k 

thermal  conductivity  of 
accelerating  gas 

Ts 

particle  sublimation  tempera¬ 
ture 
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time  it  takes  for  particle  to 
reach  its  sublimation  tempera¬ 
ture 

time  for  sublimation  to  zero 
radius,  eq.  (14) 

any  time  after  arrival  of  shock 


volume 
gas  velocity 

particle  velocity,  ft/sec 

velocity  of  particle  at  time  tj 

position  of  contact  surface 

distance  from  injection  point  of 
particle,  ft 

position  of  particle  at  time 

total  kinetic  energy  of  N  parti¬ 
cles  falling  on  area  A^, 


joules 


viscosity 

wavelength 

density  of  particle 

total  final  reflectance  of  area 
after  ejqjosure  to  e 

total  initial  reflectance  of 
area  A^  before  exposure 

total  effective  reflectance  for 
area  damaged  by  crater 

spectral  hemispheric  angular 
reflectance 

Stefan -Boltzmann  constant 
gas  flow  duration 
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APPENDIX  B 


CALCULATION  OF  TEMPERATURES  FOR  VISIBILITY  ON 
FILM  OF  HEATED  PARTICLES 

Temperature  History  of  Accelerated  Particle 

The  heat-balance  equation  between  the  gas  and  the  particle  is 

-VC  p  ^  =  Ah(T  -  T  J  (9) 

P  dt  ^ 


Over  the  temperature  range  involved,  an  appropriate  average  C  was  used.  With  the 

IT 

assumption  of  spherical  particles,  equation  (9)  becomes 


-  -  nr^Co  ^  =  4ffr^h(T  -  T  J  (Bl) 

3  P  dt  ^ 


The  experimental  data  of  reference  18  for  spheres  indicated  that  for  Re  <  10,  hD/k 
approached  1,  although  the  theoretical  value  approaches  2  for  Re  —  0.  Choice  of  a  con¬ 
stant  value  of  1. 0  for  hD/k  is  strengthened  by  the  work  of  reference  9.  Although  the  heat 
transfer  in  reference  9  was  done  for  cylinders,  the  effect  of  Mach  number  is  shown  in  the 
slip-flow  region  and  indicates  that,  as  the  Mach  number  decreases,  the  Reynolds  number 
also  decreases  in  such  a  manner  that  the  Nusselt  number  stays  approximately  constant. 
Thus,  if  h  =  k/D  is  substituted  into  equation  (Bl) 


3k 


2r2c  p 

Letting  C  ^  results  in 

^  3k 


(B2) 


(B3) 
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On  integration,  equation  (B3)  becomes 


T  H=T+(T-T)e  (B4) 

Minimum  Particle  Temperature  Required  to  Expose  Film 

The  minimum  energy  needed  to  expose  high-speed  panchromatic  film  with  a  1250  ASA 
“4 

rating  is  76x10  erg  per  square  centimeter.  The  value  was  obtained  by  the  manufacturer 
by  using  a  standard  tungsten  light  source.  Now  the  tungsten  source  radiates  at  about  the 
same  temperature  as  the  particles  in  the  hot  gas  (~3000°  K).  The  film  used  is  sensitive 
between  0.  3  to  0.  7  micron,  and  at  these  wavelengths  there  is  little  difference  in  the 
spectral  matching  to  the  film  of  radiation  given  off  by  the  particles  and  that  given  off  by 
the  calibrating  source.  Thus  the  energy  density  figure  given  was  used  unchanged.  To 
determine  the  temperature  needed  by  the  particles  to  radiate  sufficient  energy  to  expose 
the  film,  the  radiation  law  is  used  in  the  form 


EAf  =  £OA(Tp 

\47rR  / 


where  A^  is  the  area  of  the  film  exposed  by  the  particle  radiation.  Assuming  e  =  1, 
that  is,  a  blackbody,  using  the  known  demagnification  factor  (F  =  L/L^)  of  2.  5  from  the 
geometry,  and  substituting  for  A^,  A,  and  give  equation  (B5)  in  the  form 


2  /  2' 

E(7rr^)  2  4/ 

^ - i=a{Airv^){T  2)  — ^ 

(2. 5)2  V4,r2 


The  time  of  exposure  t  was  assumed  equal  to  the  time  it  takes  for  a  particle  to  pass 
through  its  diameter  D  going  at  a  velocity  V  (t); 


Vg(l  -  e"^^) 


Substituting  equation  (B7)  into  equation  (B6)  and  solving  for  (T  2)  yield 


4  ER^Vjt)  ERV(l-e"^b 
(Tn  2>  = - - - = - - - 

p,  ^  00  9  9 

(2.  5)'^arp  (2.  5)  Vp 


30 


REFERENCES 


1.  D’Aiutolo,  C.  T. :  Meteoroid  Hazards  in  Near  Earth  and  Deep  Space.  Paper  Pre¬ 

sented  at  Manned  Planetary  Mission  Tech.  Conf. ,  Lewis  Research  Center,  Cleve¬ 
land  (Ohio),  May  21-23,  1963. 

2.  Alexander,  W.  M. ;  McCracken,  C.  W. ;  Secretan,  L. ;  and  Berg,  O.  E. ;  Review  of 

Direct  Measurements  of  Interplanetary  Dust  from  Satellites  and  Probes.  Goddard 
Space  Flight  Center  Contributions  to  the  Cospar  Meeting,  May  1962.  NASA  TN 
D-1669,  1963,  pp.  39-60. 

3.  Golovin,  M.  N. ;  and  Putnam,  A.  A. :  The  Inertial  Impaction  of  Small  Particles. 

Rept.  No.  TN  58-36,  Arnold  Eng.  Dev.  Center,  June  1958. 

4.  Marble,  Frank  E. :  Dynamics  of  a  Gas  Containing  Small  Solid  Particles.  Combustion 

and  Propulsion,  R.  P.  Hagerty,  et  al. ,  eds. ,  Pergamon  Press,  1963,  pp.  175-213; 
Discussion,  pp.  213-215. 

5.  Rudinger,  George:  Dynamics  of  Gas-Particle  Mixtures  with  Finite  Particle  Volume. 

Paper  No.  65-9,  AIAA,  1965. 

6.  Mirels,  H. ;  and  Braun,  W.  H. :  Nonuniformities  in  Shock-Tube  Flow  Due  to 

Unsteady-Boundary-Layer  Action.  NACA  TN  4021,  1957. 

7.  Roshko,  Anatol;  On  Flow  Durations  in  Low-Pressure  Shock  Tubes.  Phys.  Fluids, 

vol.  3,  no.  6,  Nov.  -Dec.  1960,  pp.  835-842. 

8.  Glass,  I.  I. :  Shock  Tubes.  Pt.  I  -  Theory  and  Performance  of  Simple  Shock  Tubes. 

Rev.  No.  12,  Univ.  of  Toronto,  May  1958. 

9.  Laurence,  J.  C.;  and  Sandborn,  V.  A.:  Heat  Transfer  from  Cylinders.  Symposium 

on  Measurements  in  Unsteady  Flow,  Worchester  (Mass. ),  May  21-23,  1962,  pp. 
36-43. 

10.  Anon. :  Kodak  Plates  and  Films  for  Science  and  Industry.  Pub.  No.  9,  Eastman 

Kodak  Co. ,  1962. 

11.  Camm,  John  C. ;  and  Rose,  Peter  H. :  Electric  Shock  Tube  for  High  Velocity  Simu¬ 

lation.  Res.  Rept.  No.  136,  Avco-Everett  Res.  Lab. ,  July  1962. 

12.  Rudinger,  G. :  Experiments  on  Shock  Relaxation  in  Particle  Suspensions  in  a  Gas  and 

Preliminary  Determination  of  Particle  Drag  Coefficients.  Rept.  No.  CAL-90-P, 
Cornell  Aero.  Lab. ,  July  1963. 

13.  Crowe,  C.  T. ;  Nicholls,  J.  A.;  and  Morrison,  R.  B. :  Drag  Coefficients  of  Inert  and 

Burning  Particles  Accelerating  in  Gas  Streams.  Ninth  Symposium  (International) 
on  Combustion,  Academic  Press,  Inc. ,  1963,  pp.  395-406. 


31 


14.  Summers,  James  L. :  Investigation  of  High-Speed  Impact:  Regions  of  Impact  and 

Impact  at  Oblique  Angles.  NASA  TN  D-94,  1959. 

15.  Herrmann,  Walter;  and  Jones,  Arfon  H. :  Survey  of  Hypervelocity  Impact  Informa¬ 

tion.  Kept.  No.  99-1,  M.I.T.,  Sept.  1961. 

16.  Eichelberger,  R.  J. ;  and  Gehring,  J.  W. :  Effects  of  Meteoroid  Impacts  on  Space 

Vehicles.  BRL  Rept.  No.  1155,  Aberdeen  Proving  Ground,  Dec.  1961. 

17.  Mark,  Herman;  Sommers,  Ralph  D. ;  and  Mirtich,  Michael  J. :  Effect  of  Surface 

Thermal  Properties  of  Calibrated  Exposure  to  Micrometeoroid  Environment. 

Paper  No.  65-138,  AIAA,  1965. 

18.  McAdams,  W.  H. :  Heat  Transmission.  Third  Edition,  McGraw-Hill  Book  Co. ,  Inc. , 

1954. 


32 


NASA-Langley,  1966  E-3039 


"The  aeronautical  and  space  activHies  of  the  United  States  shall  he 
conducted  so  as  to  contribute  ...  to  the  expansion  of  human  knowl¬ 
edge  of  phenomena  in  the  atmosphere  and  space.  The  Administration 
shall  provide  for  the  widest  practicable  and  appropriate  dissemination 
of  information  concerning  its  activities  and  the  results  thereof." 

— National  Aeronautics  and  Space  Act  of  1958 


NASA  SCIENTIFIC  AND  TECHNICAL  PUBLICATIONS 


TECHNICAL  REPORTS:  Scientific  and  technical  information  considered 
important,  complete,  and  a  lasting  contribution  to  existing  knowledge. 

TECHNICAL  NOTES:  Information  less  broad  in  scope  but  nevertheless 
of  importance  as  a  contribution  to  existing  knowledge. 

TECHNICAL  MEMORANDUMS:  Information  receiving  limited  distri¬ 
bution  because  of  preliminary  data,  security  classification,  or  other  reasons. 

CONTRACTOR  REPORTS:  Technical  information  generated  in  con¬ 
nection  with  a  NASA  contract  or  grant  and  released  under  NASA  auspices. 

TECHNICAL  TRAlNSLATIONS:  Information  published  in  a  foreign 
language  considered  to  merit  NASA  distribution  in  English. 

TECHNICAL  REPRINTS:  Information  derived  from  NASA  activities 
and  initially  published  in  the  form  of  journal  articles. 

SPECIAL  PUBLICATIONS:  Information  derived  from  or  of  value  to 
NASA  activities  but  not  necessarily  reporting  the  results  of  individual 
NASA-programmed  scientific  efforts.  Publications  include  conference 
proceedings,  monographs,  data  compilations,  handbooks,  sourcebooks, 
and  special  bibliographies. 


Details  on  the  availability  of  these  publications  may  be  obtained  from: 

SCJENTIFIC  AND  TECHNICAL  INFORMATION  DIVISION 

NATIONAL  AERONAUTICS  AND  SPACE  ADMINISTRATION 

Waihinston,  D.C  20546 


